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Simulation of anyonic fractional statistics of Kitaev's toric model 
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Abstract. - Since the anyonic excitations in the Kitaev toric model are perfectly localized quasi- 
particles, it is possible to generate dynamically the ground state and the excitations of the model 
Hamiltonian to simulate the anyonic interferometry. We propose a scheme in circuit QED to sim- 
ulate the interferometry. The qubit-cavity interaction can be engineered to realize effective state 
control as well as the controlled dynamics of qubits, which are sufficient to prepare the ground 
states, create and remove the anyonic excitation, and simulate the anyonic interferometry. The 
simplicity and high fidelity of the operations used open the very promising possibility of simulating 
fractional statistics of anyons in a macroscopic material in the near future. 
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Introduction. — Anyons are exotic quasiparticles liv- 
ing in two dimensions with fractional statistics 1 . A 
paradigmatic system for the existence of anyons is a kind 
of so-called fractional quantum Hall states [2]. Alterna- 
tively, artificial spin lattice models are also promising for 
observing these exotic excitations [3H5|. Recently, with 
the potential applications in topological quantum com- 
putation, anyons have attracted strong renewed interests 
|5Ml5|. However, a direct observation of fractional statis- 
tics associated with anyon braiding is hard. For the Ki- 
taev toric model, it has been proposed an alterna- 
tive, i.e., to generate dynamically the ground state and 
the excitations of the model Hamiltonian, instead of direct 
ground-state cooling, to simulate the anyonic interferom- 
etry. The anyons are perfectly localized quasiparticles in 
that model Hamiltonian, therefore one do not need a large 
system for implementing their braiding operations. In- 
deed, small-scale system for proof-of-principle simulation 
of the anyonic fractional statistics was realized in optical 
systems 12 11. This kind of experiment represents an 
important step toward the long pursued goal to simulate 
fractional statistics of quasiparticles in a macroscopic ma- 
terial, which is critical for future implementation of topo- 
logical quantum computer, providing the ability to scale 
up to large systems. 

Superconducting circuit is one of the most promising 
candidates serving as hardware implementation of quan- 



tum computers [16j . In this paper, we propose a scheme 
in circuit QED to simulate the interferometry. The qubit- 
cavity interaction can be engineered to realize effective 
state control as well as the controlled dynamics of the 
two- level systems, which are sufficient |15j to prepare the 
ground states of the toric Hamiltonian, create and remove 
the anyonic excitation, and simulate the anyonic interfer- 
ometry. The simplicity and high fidelity of the operations 
used may open the very promising possibility of simulating 
fractional statistics of anyons in a macroscopic supercon- 
ducting phase qubit in the near future. 

Quantum dynamics in circuit QED. The super- 
conducting charge qubit considered here consists of a small 
superconducting box with excess Cooper-pair charges, 
formed by a symmetric SQUID with the capacitance Cj 
and Josephson coupling energy Ej, pierced by an external 
magnetic flux <f>. A control gate voltage V g is connected to 
the system via a gate capacitor C g . Focus on the charge 
regime, at temperatures much lower than the charging en- 
ergy [E c — 2e 2 /(C g + 2Cj)] and restricting the induced 
charge [n = C g V g /(2e)] to the range of n G [0, 1], only a 
pair of adjacent charge states on the island are relevant. 
The qubits are capacitively coupled to a transmission line 
resonator which forms a ID cavity. For simplicity, we here 
assume that the cavity has only a single mode that plays a 
role. To obtain maximum coupling strength, they are fab- 
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ricated close to the voltage antinodes of the cavity. This 
coupling is determined by the gate voltage, which contains 
both the dc contribution and a quantum part. The qubits 
are working at their optimal points, where they are im- 
mune to the charge noise and possess long decoherence 
time. In the qubit eigenbasis, neglecting fast oscillating 
terms using the rotating-wave approximation, the Hamil- 
tonian describes this scenario now takes the usual Jaynes- 
Cummings form |17j 



Hjc = uj r a^a + —a z — g (aV 



(1) 



where we have assume h = 1, v is the qubit energy split- 
ting, g is the coupling strength of a qubit to the cavity, 
uj r , a and is the frequency, annihilation and creation 
operator of the cavity field, respectively. Note that the 
qubit frequency can be tune within a large range by the 
external magnetic field. Therefore, selected qubit-cavity 
interaction can be achieved. 

The evolution operator of qubit-cavity interaction 
Hamiltonian in Eq. ([1} is given by 
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where 9 = gt, and it results in an oscillation between the 
qubit and cavity states. 

Meanwhile, driving in the form of 



e^a) e~ luldt + e* (t)ae iu3it 



(3) 



on the resonator can be obtained [17] by capacitively cou- 
pling it to a microwave source with frequency lo c i and am- 
plitude e(t). Depending on the frequency, phase, and am- 
plitude of the drive, different logical operations for qubit 
can be realized. 

To get fast gate, we work with large amplitude driving 
fields, where quantum fluctuations are very small compare 
with the drive amplitude, and thus the drive can be con- 
sidered as a classical field. In this case, it is convenient 
to displace the field operators using the time-dependent 
displacement operator [15] : 

D(a) — exp (aa) — a*a) . 

Choosing ia = ui r a + e(t)e~ lWdt to eliminate the direct 
drive on the resonator, by h, from the effective Hamil- 
tonian, then the displaced Hamiltonian for a qubit reads 

El 

H D = uj r a)a + ^a z - g [(a + a) a + + H.c] . (4) 

When the drive amplitude is independent of time, and 
change to a frame rotating at a><£, the displaced Hamilto- 
nian reads 



r t A A / 

H KF = da , a+ —a z + —a x - g { 



aa + + o'er ) 



(5) 



where 6 = U) r — Ud, A = v — ujd r , and SI = 2ge/5 is 
the Rabi frequency. In the dispersive regime S 3> g, the 
effective Hamiltonian is [17] 



H x = 8a) a + 



A + g 2 /S. 



CI. 



(6) 



By choosing A + g 2 /S — 0, the Hamiltonian ^ evolves 
as a rotation around the x axis. The gate speed scales as 
t x ~ which have already been experimentally verified 

m 

We now turn to consider that the drive is sufficiently 
detuned from the qubit |A| ^> ft. Then, the effective 
Hamiltonian from Eq. §5§ is [T7] 



6a*a+ ^a z 



(7) 



which generates rotations around z axis at a rate x = 
A+g 2 /<5+f2 2 /(2A), and the gate speed scales as t z ~ 1/x- 

From the above analysis, we can see that switch on and 
off the qubit-cavity coupling can be achieved by tuning 
the external driven field. In this way, we can obtain indi- 
vidual manipulation of qubits. It worth noticing that the 
dispersive regime in single-qubit operation is induced by 
the detuning with respect to the driving, which will not 
slow the gate speed as in the dispersive coupling of qubits 
via virtual excitation of the cavity mode |17] . Specifically, 
t x ~ 1/g when e ~ 5/2. 

Note that the iSWAP gate between the cavity and jth 
qubit can be obtained from evolution operator ^ when 
t = it/ (Jig). Together with single-qubit rotations, it can 
be used to generate cluster states more efficiently than 
that of controlled phase flip gate 21) 2J]. The operation 
used for the cavity and qubit system is 



U c 
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where is a ir/2 rotation around the z axis of the 

qubit state. This gate makes the two involved parties 
linked by a controlled phase flip gate and at the same time 
it swaps their states. It is worth noticing that both single- 
qubit and iSWAP operation times are several ns which is 
much smaller compare to typical decoherence time of the 
qubit [22] and cavity [23] , both of them being on the order 
of /is. 

Ground states preparation. — In the following, we 
show how to prepare the ground states of the toric model 
[4] with this gate. The toric code is defined as the ground 
level of a stabilizer Hamiltonian 



H 



(9) 
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Fig. 1: (Color online) The toric model, (a). It is denned as the 
ground state of a stabilizer Hamiltonian on a square lattice with 
spins at the edges of a square lattice. The commuting vertices 
(plaquettes), as indicated by a blue cross (red rectangle), are 
the stabilizers. A single a x (a z ) gate on a spin, indicated by a 
small red (blue) circle, can create a pair of magnetic (electric) 
defects excitation, indicated by a big red (blue) circle, (b). 
The smallest system for implementation of the anyon braiding 
operation, the ground state of which is equivalent under local 
single-bit operations to a graph state in (c). 




Fig. 2: (Color online) The process of generating the ground 
state graph with six spins. Each dashed circle denotes a U c 
operation and each bond means the two are linked by a con- 
trolled phase flip gate. Circle with capital letter "C" and num- 
bers denote the cavity and spins, respectively. The blue arrow 
across the circle represents "X" measurement on the cavity. 



on a square lattice with spins, realized as superconducting 
qubits here, at the edges of a square lattice as shown in Fig. 
(jTJi) . The sum is over the mutually commuting stabilizers 

A v = Y[af, B p = l[a*, 

where v runs over all vertices and p over plaquettes, as 
indicated by a blue cross and a red rectangle in Fig. (QJi). 
The ground state \<p) g is characterized by A v = B p = 1. 
Plaquette and vertex excitations \(p) etTn , also called mag- 
netic and electric defects and indicated in Fig. ([Th.) by 
bigger red and blue fulled circles, are characterised by 
Bp = —1, A v = — 1 and they are mutual Abelian anyons. 
Therefore, apply single a x or a z can create these quasipar- 
ticle excitations. These excitations always appear in pairs, 
at the ends of strings of a x and a z operators applied on 
a ground state. To apply a x or a z to a spin, we need to 
turn on the qubit-cavity interaction in Eq. © or (|7|) . 

For the toric model, the excitations are perfectly lo- 
calized. Therefore, for a proof-of-principle simulation of 
anyons and their braiding statistics, the minimum imple- 
mentation only needs six spins jllj , as the graph shown in 
Fig. QJb) . The ground state of which is equivalent under 
local single-bit operations to a graph state in Fig. [T{c). 
We next show how to prepare the ground state. We work 
with a appropriate 3-body vertex and plaquette operators 
along the boundary providing for a two-dimensional code 
space. We start with a well defined state, i.e., the cavity 
and spins 1, 2, 3, 5, and 6 are initialized to the state of 
1+) = (|0) + |l))/v / 2 while spin 4 is prepared in the ground 
state \g). The process is illustrated in Fig. ([2]). We first 
sequentially apply U c on cavity and spins 2, 1, 3, 6, and 
5, which results in the five qubit cluster state and leave 
the cavity to the |0) state. As illustrated in the Fig. ([2]), 
after each U c , the cavity is at the right end of the cluster. 



Then, apply U c on spin 4 to create an bipartite entangle- 
ment. At last apply U c on spin 6 to fuse the two entangled 
states followed by an X measurement on the cavity state. 
After these steps, the six spins are now prepared in the 
graph state as in Fig. ((Hi). Measurement on the cavity 
can be implemented by swapping the state of the cavity 
to another spin and then measure the spin state. 

Controlled dynamics and anyon interference. 

This operation in Eq. can also be used to realize con- 
trolled rotation U z = |0)(0| ® a z + \1){1\ ® I of the qubit 
about the z axis under the control of the cavity states. It 
goes as following: 1) Prepare the cavity state to the state 
of (|0) + |1)). Here, we need to engineer the cavity num- 
ber states, which cab be achieved by swapping an ancillary 
qubit states with that of the cavity. 2) Excite the qubit 
excited state |e) to an ancillary level (other than \g) and 
|e)). Ancillary levels have already been used in a recent 
experimental demonstration of gates in circuit QED [22 . 
3) Apply U c for the cavity and the target spin for a time 
ofn/S. Note that U x = |0)(0|®cr a .+|l)(l|® J is equivalent 
to U z up to local single-qubit operation on the spin [10] : 
U x = HU Z H . It is note that our controlled operations are 
different from that of Refs. [10,15 . The essential of the 
controlled operations is to obtain state-dependent dynam- 
ics, in this sense, the effect of our controlled operations are 
the same as their's. 

With these controlled operators and rotating gates of 
the qubits, it is sufficient to create and remove the ex- 
citations and to simulate the anyonic interferometry of 
the toric model. To create a superposition of the ground 
and excite state, we need an ancillary degree of freedom, 
realized by the cavity here, then conditional excitation 
~ (|y) 9 |l) +^|y)e,m|0)) can be achieved, where r\ is deter- 
mined by the relative amplitude of the controlled gate. It 
is obvious that r\ = 1 is the controlled operation U X)Z . To 
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Fig. 3: (Color online) The minimal anyonic interferometry. 
The label of spins (yellow small circles) is indicated in (a). 
Applying a a z gate to spin 3, creates two electric defects e\ 
and 62, indicated as big blue circles in (b). With the cavity 
in state |+), applying a U x gate on spin 4 conditionally create 
two magnetic defects mi and m,2, i.e. the system now in the 
superposition of (c) and (d). With the cavity in state |+), ap- 
plying U x gates on cavities and spins 6, 5, 3 and 4 sequentially 
result in the superposition of (c) and (e), which is equivalent 
to (f). 



apply U z to a spin, we can get controlled magnetic defects 
excitation: U x \ip) g \+) ~ (|^) 9 |1) + \(f) m \0)). Similarly, 
one can also create controlled electric defects excitation 
by U z \<p) g \+) ~ (\<p) g \l) + \<p) e \0)). 

Once the toric code state is prepared, one can simulate 
the fractional statistical phase of the anyons through a 
Ramsey- type interference experiment. The simplest any- 
onic interferometry [TT] is shown in Fig. Q. Apply a a z 
gate to spin 3, creates two electric defects e\ and ei. With 
the cavity in state |+), apply a U x gate on spin 4 creates 
the superposition states of (c) and (d) . Apply U x gates on 
cavities and spins 6, 5, 3 and 4 sequentially will move mi 
around ei and finally fusing with ui\ to vacuum. By this 
braiding, the -1 factor appears in the cavity state, i.e. the 
system now in the superposition of (c) and (e), which is 
equivalent to (f). In this interferometry the cavity state 
will change to |— ) = H\l) from the state of |+). This 
phase factor on the cavity state is solely due to the mu- 
tually fractional statistics of anyons and can be detected 
unambiguously in experiment. This detection method is 
similar with that of Ref . [TS] . 

Conclusion. — In summary, we propose a scheme in 
circuit QED to simulate the fractional statistics of anyons 
of Kitaev's toric code model. The qubit-cavity interac- 
tion can be engineered to realize effective state control as 
well as the controlled dynamics of qubits, which are suf- 
ficient to prepare the ground states, create and remove 
the anyonic excitation, and simulate the anyonic interfer- 
ometry. The simplicity and high fidelity of the operations 
used here open the very promising possibility of simulating 
fractional statistics of anyons in a macroscopic material in 
the near future. 



I thank Prof. Shi-Liang Zhu for many helpful sug- 
gestions. This work was supported by the NSFC 
(No. 11004065), the NSF of Guangdong Province (No. 
10451063101006312), and the Startup Foundation of 
SCNU (No. S53005). 

REFERENCES 

[I] F. Wilczek, Phys. Rev. Lett. 48, 1144 (1982). 

[2] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das 

Sarma, Rev. Mod. Phys. 80, 1083 (2008). 
[3] X.-G. Wen, Phys. Rev. Lett. 90, 016803 (2003); Quantum 

Field Theory of Many-Body Systems (Oxford University 

Press, Oxford) 2004. 
[4] A. Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003); ibid. 321, 2 

(2006). 

[5] Z.-Y. Xue, S.-L. Zhu, J. Q. You, and Z. D. Wang, Phys. 

Rev. A 79, 040303(R) (2009); Y. Yu and Y. Li, J. Phys. 

A: Math. Theor. 43, 105306 (2010). 
[6] L.-M. Duan, E. Demler, and M. D. Lukin, Phys. Rev. Lett. 

91, 090402 (2003). 
[7] A. Micheli, G. K. Brennen, and P. Zoller, Nature Phys. 2, 

341 (2006). 

[8] J. Q. You, X.-F. Shi, X. Hu, and F Nori, Phys. Rev. B 81, 
014505 (2010). 

[9] C. Zhang, V. W. Scarola, S. Tewari, and S. Das Sarma, 
Proc. Natl. Acad. Sci. U.S.A. 104, 18415 (2007). 

[10] L. Jiang, G. K. Brennen, A. V. Gorshkov, K. Hammerer, 
M. Hafezi, E. Demler, M. D. Lukin, and P. Zoller, Nature 
Phys. 4, 482 (2008). 

[II] Y.-J. Han, R. Raussendorf, and L.-M. Duan, Phys. Rev. 
Lett. 98, 150404 (2007). 

[12] C.-Y. Lu, W.-B. Gao, Otfried Giihne, X.-Q. Zhou, Z.-B. 

Chen, and J.-W. Pan, Phys. Rev. Lett. 102, 030502 (2009). 
[13] J. K. Pachos, W. Wieczorek, C. Schmid, N. Kiesel, R. 

Pohlner, and H. Weinfurter, New J. Phys. 11, 083010 

(2009). 

[14] J. Zhang, C. Xie, K. Peng, and P. van Loock, Phys. Rev. 

A 78, 052121 (2008). 
[15] M. Aguado, G. K. Brennen, F. Verstraete, and J. I. Cirac, 

Phys. Rev. Lett. 101, 260501 (2008). 
[16] J. Q. You and F. Nori, Phys. Today 58, (11) 42 (2005). 
[17] A. Blais, J. Gambetta, A. Waflraff, D. I. Schuster, S. M. 

Girvin, M. H. Devoret, and R. J. Schoelkopf, Phys. Rev. A 

75, 032329 (2007). 
[18] M. O. Scully and M. S. Zubairy, Quantum optics (Cam- 
bridge University Press, New York, 1997), chapt. 2. 
[19] A. Wallraff et al., Nature (London) 431, 162 (2004). 
[20] T. Tanamoto, Y.-X. Liu, X. Hu, and F. Nori, Phys. Rev. 

Lett. 102, 100501 (2009). 
[21] L.-B. Yu, Z.-Y. Xue, Z. D. Wang, Y. Yu, and S.-L. Zhu, 

larXiv:0904.T275l (2009), to be published in Eur. J. Phys. 

D, DOI: 10.1140/epjd/e2010-00258-5. 
[22] L. DiCarlo et al., Nature (London) 460, 240 (2009). 
[23] L. Frunzio, A. Wallraff, D. I. Schuster, J. Majer, and R. J. 

Schoelkopf, IEEE Trans. Appl. Superconductivity 15, 860 

(2005). 



p-4 



